N = [2:65,70:5:200]';
x = load('tableData.dat');

Hnumerical = x(:,2) ./ 2;

o = optimset('MaxFunEvals',4000,'MaxIter',4000,'Display','off');

x1 = lsqcurvefit(@fitFunctionForHWithLeadingOrder,[0,0,0,0,0,0,0],N,Hnumerical,[],[],o);

x1